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ABSTRACT 

We have used the CANGAROO-III imaging atmospheric Cherenkov tele- 
scopes to observe the high-frequency-peaked BL Lacertae (HBL) object 
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PKS 2155-304 between 2006 July 28 (MJD 53944) and August 2, triggered 
by the H.E.S.S. report that the source was in a high state of TeV gamma-ray 
emission. A signal was detected at the 4.8 a level in an effective live time of 
25.1 hours during the outburst period. The flux of Very High Energy gamma- 
rays from the CANGAROO-III observations shows the variability on the time 
scale of less than a few hours. The averaged integral flux above 660 GeV is 
(1.6 ± 0.3stat ± Q.^syst) X 10~^^ cm~^ scc"^ which corresponds to ~ 45% of the 
flux observed from the Crab nebula. Follow-up observations between August 17 
(MJD 53964) and 25 indicate the source activity had decreased. 

Subject headings: Galaxies: active — BL Lacertae objects: Individual: PKS 2155— 304 
— Gamma rays: observations 



Introduction 



The most astonishing contribution from the ground-based lACTs has been the detec- 
tion of AGNs at energies above several hundred GeV. To date, 19 active galactic nuclei 



have been reported to emit Very High Energy (VHE) gamm a-ravs (iHinton 



2007 



most 



of them classified as high-frequency-peaked BL Lacs (HBLs) (jPadovani fc Giommil Il995l ). 
Some, such as Mrk 421 and Mrk 501, have been targets of simultaneous n i ulti- w avelength 
campaigns since their first TeV detection ( Punch et al. 19921 : Quinn et al. 1996), and are 



well studied. The derived broadband spectral energy distributions (SEDs) show that they 
have two components. The lower energy component is believed to be synchrotron emission 
produced by relativistic electrons in the jet, and the higher energy component, although 
relatively poorly understood, is probably explained by the inverse- Compton scattering of 



(see e.g.. 


Jones. O'dell & Stein 


1974; 


Ulrich. Maraschi. & Urrv 


1997 


Dermer & Schlickheiser 


1993; 


Sikora. Beeelman & Rees 


1 199^ 


D- 



Multi-wavelength observations of HBLs have also shown that the fluxes are extremely 
variable, particularly in the higher energy bands. Flux variations on a wide range of time 



scales from months to two minutes are reported at TeV energies fe.g.. iGaidos et al 



Krawczvnski et al. 2001 



1999 



Gliozzi et al. 



2006 



Aharonian et al.ll2002l . l2005d : iBlazejowski et al.ll2005l : lAharonian et al 



1996 



Albert et al.ll2007al ). Time variability is a very useful way to probe 
the emission mechanism in the jet. Characteristic time scales, which can be employed as 
a measure of super-massive black hole mass at the center of HBL, have been studied with 
X-ray data using stati s tical techniques su c h as power spectra, structure fu nctions, etc (e.g., 
Havashida et all Il998l : Izhang et all Il999l : iKataoka et aD I2OO1I : ICui I I2OO4I ). A good corre- 
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lation of the X-ray and the TeV fluxes has also been reported (e.g 



Takahashi et al 



1996 



Krawczynski et al 



Diannati-Atai et al.lll999: Takahashi et al. 2000 



2002 



. Buckley et al. 


1996; 


Sambruna et al. 


2000 



Albert et al.ll2007al ). Recently, however, poor correlations or a total 



lack of correlation between the X-ray and TeV fluxes o f individual flares have been observed 



[e.g., IBlazejowski et al.l 120051 : lAtiaroman et al.l 12005a). and even TeV flares wit h no X-ray 
counterpart have been found (IBlazejowski et al.l l2005l : iKrawczynski et al.ll2004r). Spectral 



hardening with flux in the TeV energy range has been reported f or Mrk 421 (lAharonian et al. 



2002; Krennrich et al. 2002; Albert et al. 



2007a|) and Mrk 501 flAlbert et al.ll2007bh . 



X-ray (IBrinkmann et al. 
I999I 



1994 



Kubo et al. 



P KS 2155-304 (z = 0.116 flFalomo et al.lll993r)1 is one of the brightest BL Lacs in the 



199 



Giommi et al 



1998 



Vestrand fc Sreekumar 



Nicastrol I2OO2I : Izhang et all boosh and EUV jMarshall. Carone. fc Fruscionel Il993h 



bands. Since the discovery of X-ray emission from this object (iGriffiths et al. I I1979I ; ISchwartz et al.l 



19791 ). it has been repeatedly observed o ver a wide range of 



High 



1995 



nergy (VHE) gamma-rays (e.g., Treves et al. _ 19891; 



Zhang fc Xie 



Aharonian et al.ll2005al ) 



1996; Plan et al. 1997; Finer fc Edwards 



requencies from radio to Very 



Vestrand. Stacy, fc Sreekumar 



2004 



Dominici et al 



2004 



In 1997 Noven iber, a gamma-ray and X-ra y outburst fr om FKS 2155— 304 wa s de- 
tected by EGRET JSreekumar fc Vestrandl ll997h. BeppoSAX Jchiappetti et al.lll999h and 
RXTE (jVestrand fc Sreekumarlll999l ). During this active phase, the Durham group observed 
FKS 2155—304 using their Mark 6 tel escope, and reported th e flrst detection of VHE gamma- 
rays at the 6.8 a level above 300 GeV (IChadwick et al. I ll999ah at the time of BeppoSAX obser- 
vatior i. However they found no evidence for TeV gamma-ray emission during observations in 
1998 (IChadwick et al.lll999bl ) when the X-ray flux level was low. FKS 2155—304 was observed 
with the CANGAROO-I 3.8 m telescope in 1997, though poor weather resulted in minimal 
overlap with the Du rham observations. No gamma-ray signal above 1.5 TeV was detected 
jRoberts et al.lll999h . FKS 2155-304 was further observed in 1999, 2000 and 2001 with the 
CANGAROO-II telescope. It remained in a low state of X- ray activity in those p e riods, and 
was not detected above the energy threshold of 420 GeV (INishijima et al.l I2OOII ; iNishijima 
2OO2I ; iNakase et al.ll2003l ). FKS 2155—304 was conflrmed as a TeV gamma-ray source by the 
H.E.S.S. group in 2004. They reported a 45 a detection at energies greater tha n 160 GeV 
in July and October, 2002, and June-September, 2003 (lAharonian et al.ll2005al ). The flux 
variability on time scales of months, days, and hours were also reported, and the monthly- 
averaged flux above 300 GeV was between 10% and 60% of Crab flux. The energy spectrum 
was characterized by a steep power law with a time-averaged photon index of F ~ 3.3. Multi- 
wavelength observations in 2003 in a low state showed no correlation between the X-ray and 
the gamma-ray fluxes, or between any other wavebands, eve n though the fastest eve r X-ray 
flare in this object, with a 1500 sec timescale, was detected (lAharonian et al.ll2005bl ). 
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In July 2006, the H.E.S.S. group reported that PKS 2155-304 had been detected a t 
historically high TeV flux levels of up to several crab (IBenbow et al.ll2006l : iRaue et al.ll2006l ). 
They also presented very rapid flux variability with one-minute ti me scale resolution an d 
found well-resolved bursts varying on timescales of ~ 200 seconds (lAharonian et al.l 120071 ). 
The H.E.S.S. report triggered multi-wavelength Target of Opportunity (ToO) observations, 
including the observations with the CANGAROO-III telescopes from 2006 July 28 (MJD 
53944) described in this paper. 

Numerous mul t i-wavelength campaigns of PKS 2155— 30 4 have been carried out (e.g.. 



Edelson et al.lll995l : lUrrv et al.lll997 



tions for the emission mechanisms have been reported fe.g..lChiappetti et al 



2000l : lEdelson et all 1200 ll : iTanihata et al.l l200ll : IZhang et al. 



Aharonian et al.ll2005bl'). and many physical implica- 



200^ 



i l2006h . 



19991 : iKataoka et al. 
However, these 



models are still face challenges in explaining the complex and different patterns of multi- 
wavelength variations for each epoch. 

For objects like PKS 2155— 304 which show rapid and complex time variability, contin- 
uous monitoring is very important for modeling the emission mechanisms. The difference 
in longitude between the H.E.S.S. and CANGAROO-III sites is ~ 120°, corresponding to 
an 8 hour time difference. Thus the H.E.S.S. and the CANGAROO-III data complement 
one another. Northern hemisphere blazars, such as Mrk 421, Mrk 501 and lES 1959-1-650, 
have been observed at T eV energies continuously with more thari two geographically distan t 



telescope systems (e.g., iRebillot et al.l l2006l : iGliozzi et al.l l2006l : iKrawczynski et al.l 120041 ) 



providing greatly improved time coverage, but PKS 2155—304 is the first object in the 
southern sky for which such studies have been made. 



Detector and Observations 



The CANGAROO-III imaging atmospheric Cherenkov telescope system which consists 
of four telescopes are operated in Woomera, South Australia (longitude 136°47'E, latitude 
31°06'S, 160 m a.s.l.). They have parabolic reflectors of 10 m diameter with an 8 m focal 
length. Each reflector consists of 114 small spherical mirrors of 80 cm diameter, which are 
made of Fiber Reinforced Plastic. The imaging cameras have 427 pixels of 0.17° size and 
a total field of view of ~ 4°. These telescopes use altitude-azimuth mounts and are placed 
at the corners of a diamond with 100 m sides. The details of the mir ror, the came r a, a nd 



the total performaii ce of the light coUecti r ig syst em are described in iKawachi et al.l (120011 ) 
Kabuki et aP J2003h . and lEnomoto et all J2006al ). 



The CANGAROO-III observations of PKS 2155— 304 during this outburst period were 
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made on five moonless nights between July 28 (MJD 53944) to August 2 in 2006. There were 
no observations on August 1, and the July 29 observations were affected by cloud. Three 
(T2, T3, T4) of the four telescopes were used in these observations. For each telescope, 
when the number of hit pixels exceeds four within a 20 nsec coincidence window, a local 
trigger is sent to the stereo trigger system through an optical link. In the stereo trigger 
system, determination of stereo events is done by requiring that at least two local trigger 
signals coinciding for at least 10 nsec within a 650 nsec time window, taking the geometrical 



time delay, which depends o n the telescope pointing direction, into consideration (IKubo et al. 



200ll : iNishiiima et al.ll2005r ). 



These observations were made using wobble mode in which pointing position of each 
telescope was shifted in declination between ±0.5° from the center of PKS 2155— 304 alter- 
natively every twenty minutes. Due to a mechanical tracking problem with the T3 telescope 
during this period, stereo observations with three telescopes were done only after culmina- 
tion. Stereoscopic observations using T2 and T4, which are located on the long diagonal of 
the diamond, were performed before culmination from July 28 to July 31. 

The details of the CANGAROO-III observations of PKS 2155-304 are summarized in 
Table [H The typical trigger rates, the mean zenith angle, and the total observation time 
of two-fold coincidences of T2-T4 were ~20 Hz, 26.6°, and 11.4 hours, respectively. The 
typical trigger rate of three-fold coincidences was ~12 Hz, with a mean zenith angle for 
these observations of 20.4° and a total observation time of 17.6 hours. The daily observation 
time tabs, average zenith angle z, and trigger rate r^^ are shown in Tabled] for both T2-T4 
and three-fold data. The zenith angles of observations ranged from 52° to 7° for T2-T4 data, 
and from 1° to 46° for three-fold data. 

Follow-up observations were made with the same system on six moonless nights between 
August 17 (MJD 53964) and 25. These data were taken only after culmination using three- 
fold coincidences. The zenith angle during these observations ranged from 1° to 48° with 
a mean of 20.9°, and a total observation time of 19.1 hours. These observations are also 
summarized in Table [TJ 



3. Analysis 

In order to calibrate the relative gain and the timing for each pixel, the data taken 
by illuminating the focal plane uniformly with LED photons were used. Conversion factors 
from ADC value to photoelectrons for each pixel are determined using the LED system in 
the camera vessel, and the timing jitter of the pulse signals are calibrated using light from a 
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LED mounted at the center of t he telescope di s h. De tails of the standard CANGAROO-III 
calibration methods are given in iKabuki et al. J2003h . These LED runs were performed at 
least once per night. Five out of 427 pixels for T2 and two for T3 were identified as bad 
pixels which gave no signal or false signals, and these were removed from the analysis. There 
were no bad pixels in the T4 camera. 

We have cleaned the images to eliminate night sky background photons and select 
shower images. Each pixel is required to have a signal larger than 5.0 p. e. with an arrival 
time within ±30nsec of the average shower arrival time. Images containing at least five 
adjacent pixels which pass the above requirements are selected. The distributions of shower 
rate each minute are shown in Fig. [1] for T2-T4 coincidence (left panel) and three-fold 
coincidence (right panel), respectively, and the average shower rate r^h for each night is listed 
in Table [H In order to ensure reliable arrival direction and energy estimations, the data 
taken during periods when shower rate was lower than ~5 Hz for three-fold coincidences and 
~10 Hz for T2-T4 coincidences are not used, as they are probably affected by clouds. In 
the following analysis, we further require that none of the brightest 15 pixels of each image 
should be in the outermost layer of the camera in order to avoid deformation of the image. 
This met hod improves the effect ive area by more than 20% compared to the previous simple 
edge cut (lEnomoto et al.ll2006al ). particularly above ~2 TeV. After taking into account the 
DAQ dead-time, the effective total live time tuv is also summarized in the same table. 

The monients of the shower image are then parameterized using the Hillas parameters 
( iHillas et al.lll985l ). and the arrival directions are reconstructed using the intersection of 
image axes. The intersection point is obtained by minimizing the sum of squared widths of 
the images seen from the assumed point with a constraint on the distances between images' 
cent er of gravity a n d assu med intersection point (IP-fit) considering the len gth I width ratio 



(see iKabuki et al.l (120071 ) for details), which is similar to Algorithm 5 in iHofmann et al. 
f|l999f )~ This method i mproves the signal to n oise ratio by more than 10% compared to the 



procedure described in lEnomoto et al.l (l2006a] ) and is confirmed by Monte Carlo simulations 
and observations of the Crab nebula. 

After event reconstruction, numero us cosmic ra y background events are rejected using 
the Fisher Discriminant (FD) method (lFisherl ll936l) . The FD value for each event is cal- 
culated as explained in lEnomoto et al. J2006bl ) using the size-corrected width and length 
for each telescope, where the energy dependence of the width and length are corrected to 
keep each parameter at the same values in various energies. In order to determine the opti- 
mum FD cut values, a Monte Carlo simulation was carried out assuming a single power-law 
spectrum, dF/dE oc E~^, with an index of F = 3.3, as determined for PKS 2155—304 by 



the H.E.S.S. group (lAharonian et al.l l2005al ). Before applying the FD cuts, the data are 



- 7- 



divided into two data sets, " 2 fold" which include the data taken by only T2-T4 coincidence 
and "3 fold" which required a three-fold coincidence between the three telescopes. Each 
dataset is further divided by zenith angle into two datasets, "Sfold-SZ" and "2fold-SZ" 
for z < 30°, and "3fold-LZ" and "2fold-LZ" for z > 30°, respectively. Considering the 
Point Spread Function (PSF) of 0.23° (FWHM), the 6^ (where 6 is the angular difference 
between the reconstructed arrival direction and the source position) cuts were applied at 
6*^ < 0.06 deg^ for the 3fold data, and to keep the same efficiency for gamma rays as for the 
three-fold analysis, 9'^ < 0.14 deg^ was applied to the 2 fold data. Based on the Monte Carlo 
simulation assuming 100 % Crab flux of gamma rays above 600 GeV, we could determine 
the best FD cut criteria which yielded the highest significance for each dataset. Following 
the results from the Monte Carlo study, we select the events with FD > —0.3 for 3fold- 
SZ, FD > -0.2 for 3fold-LZ, FD > -0.1 for 2fold-SZ and FD > -0.2 for 2fold-LZ 
datasets, as candidate gamma-ray events. 

The primary gamma-ray energy is estimated from the number of detected photo-electrons, 
based on the Monte Carlo simulations assuming a single power-law spectrum as mentioned 
above. This relation depends on the zenith angle of observations, so the simulations were 
done using the same variation of elevation as the actual observations. Although there is also 
some dependence on the impact parameter in the energy determination, we do not incorpo- 
rate that here, with the resulting energy resolution estimated to be 30 % around 1 TeV. The 
effective detection area has been estimated from the Monte Carlo simulation as a function 
of energy. The detection energy threshold Efh is taken to be the energy of the peak of the 
distribution of triggered shower energies. 

4. Results 

In this section, we show the results of CANGAROO-III observations of PKS 2155— 304 
in 2006. A summary of the results are shown in Table [2j Using these results, we present the 
average integral flux, the average differential energy spectrum, and the time variation. 

4.1. Average flux 

After the data reduction described in previous section, the flnal 9"^ distribution is ob- 
tained. Fig. [2] (a) and (c) show the results in the outburst period for 3fold and 2fold 
datasets, respectively. The background level was estimated using off-source data in the cor- 
responding region in the other side of the fleld of view. The distributions of the excess events 
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against i.e., with the off-source events subtracted from the on-source events, are shown 
in Fig. Et^b) and (d) for "if old and 2fold, respectively. There are clear excesses at small 
values of 6"^ in both figures corresponding to the observed signal from PKS 2155— 304. The 
number of excess events are calculated using a circular region of 9^ < 0.06 deg^ for 3 fold and 
6*^ < 0.14 deg^ for 2 fold centered on the source position, respectively, considering the PSF as 
mentioned in the previous section. The hatched histograms in figure [2] (b) and (d) indicate 
the expected PSF normalized to the number of excess events in 6^ < 0.25. The combined 
excess of all the data is 322±67 events (4. 80"), where only the propagation of statistical errors 
is considered. Fig. [3] shows the FD distributions after the 6^ cut for on-source events, off- 
source events and the gamma-ray candidate events (on-source— off-source) together with the 
Monte-Carlo gamma-rays normalized to the number of excess events. The FD distributions 
are quite consistent for the observed gamma-ray events and the Monte Carlo gamma rays. 
The FD cut is not applied for this plot. The time- averaged integral flux above 660 GeV is 
calculated to be F{> 660 GeV) = (1.6 ± 0.3,^* ± 0.5syst) x 10"^^ cm"^ sec"^ Thi s cor- 



responds to ~ 45% of the flux observed from the Crab nebula (lAharonian et al.ll2004r). and 



is a fa ctor of five more intense than the flux in the low state reported in (lAharonian et al. 



2005al ). The systematic error in the flux arises mainly from energy scale uncertainties due 
to the absolute light collection efficiency (20%), FD cut criteria (22%), uncertainties in the 
probability density function of images generated by Monte Carlo simulation (10%), and the 
uncertainty of a power law index (8%). At a higher energy threshold above 1.0 TeV, which 
corresponds to the energy threshold for 3fold-LZ and 2fold-LZ data, no significant excess 
above our sensitivity from PKS 2155— 304 is found in this period, and the resulting 2a flux 
upper limit is F{> 1.0 TeV) < 9.2 x 10"^^ cm'^ sec'^ 

For a cross-check, an alternative analysis using the Fisher Discriminant with a fit 



( lEnomoto et al.l l2006bl ) was performed on the same data. Although this is our standard 
method of analysis, the fine time binning and consequent small count statistics make it 
difficult to estimate the systematic errors correctly. The main source of error is the non- 
uniformity of the Fisher Discriminant in the ring-shaped background region. So we primarily 
took a conservative FD cuts method using a single reflected background at the expense of 
the statistics, and used FD fits as a cross-check. The difference of the integral flux between 
them was less than 15% comparing at the same threshold energy. 



4.2. Average differential energy spectrum 



We obtained the time-averaged differential energy spectrum from all the data except for 
large zenith angle dataset, i.e. from 3fold-SZ and 2fold-SZ, since the large zenith angle 
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data have a higher energy threshold and therefore lower energy bins have a different exposure 
than higher energy one. The differential flux is shown in Fig. H] as closed circles with la 
statistical error bars. The best fit of a power law to the small zenith angle data yields a 
photon index F = 2.5 ± O.bgtat =t O.Tgyst and a flux normalization A^o( 1 TeV) = (1.0 ± 
0.2stat ± O.Ssyst) X 10"^^ cm~^ sec~^TeV~^. We estimate the systematic error considering the 
same factors as mentioned in the previous section. 



4.3. Time variations 

We have searched for gamma-ray emission from PKS 2155— 304 on a night-by-night 
basis. The hve time, tuy, the number of on-source and off-source events, Nqn and Nqff, the 
number of excess events, A^, and significance, s, corresponding to each night are summarized 
in Table [2l Here the off-source data are summed up over the five nights to reduce the 
statistical fluctuation, and normalized to the live time of each night. The nightly average 
integral fluxes F above 660 GeV, using all the data are also shown in the same table, and 
are plotted in Fig. [5l where the flux from 3 fold and 2 fold datasets are combined. The light 
curve shows that the average flux reached ~ 70 % Crab in the night of July 30. Assuming 
a constant average flux, a fit yields a value of 13.9 for 4 degrees of freedom, which 
corresponds to a probability of ~ 0.8 %. This means that there is a marginal time 
variations of the average flux on a nightly basis between July 28 and August 2. 

We further divide the data into ~ 40 minutes bins for each night, and intra-night 
variation is investigated. Since the pointing offset from the source position is alternated in 
sign every 20 minutes in our observations, 40 minutes bin is adopted as it offers a reasonable 
cancellation of any asymmetric effects. The light curve of PKS 2155— 304 expressed by the 
integral flux of VHE gamma-rays above 660GeV using 2foldSZ and 3fold-SZ dataset is 
shown in Fig. [6l 

Assuming a constant average flux for the data above 660 GeV, a fit yields a value of 
96.4 for 30 degrees of freedom, which corresponds to a probability of 7 x 10"^. The same 
calculation was done for the data of July 28 and 30, and gave values of 29.2 and 22.1 for 
6 degrees of freedom, respectively. The intra-night variations are apparent. From figure [6l 
altho ugh it is difficult to calculate the fractional root mean square variability amplitude 



Fyar (jVaughan et al.l 120031 ) because of the poor statistics, Fyar = 0.75 ± 0.07 on July 28 and 



0.58 ± 0.08 on July 30 are obtained. 

The results in the follow-up observations between August 17 and 25 are also summarized 
in Tabled The 2a upper limit above 680 GeV is F{> 680GeV) < 6.8 x lO^^^ ^^-2 sec-i(< 
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20% Crab), which shows TeV gamma-ray activity had subsided compared to the outburst 
period. 



Discussion 



From Figure the fluxes on July 28 and 30 exceed 50% of the flux from the Crab 
nebula, and are the same level a s that at the end of J uly 27-28 (MJD 53944) observation 
reported by the H.E.S.S. group (lAharonian et al.l 120071 ). taking into account the difference 
of energy threshold and assuming a power law spectrum. The flux in July 29 might be 
underestimated due to the presence of clouds. Therefore, our light curve suggests that the 
flux continuously decreases on average following the large fla. r e of J uly 28. The preliminary 
H.E.S.S. light curve above 200 GeV presented by iRaue et al.l (120061 ) indicates that a 5 crab 
flare occurred between our observations of July 29 and 30. However we did not detect a flare 
in this period, possibly due to the 10 hour gap between the end of the H.E.S.S. observations 
and the start of our observations. For the shorter variability time scales, the July 28 and 30 
data indicate the intra-night variations. 

The derived values of the fractional root mean square variability amplitude, F^ar = 
0.49 for night by night variation and F^jar = 0.47 for intra-night variation in July 28, are 
comparab le to the values from the intra-night variability earlier on July 28 reported by 
H.E.S.S. ( lAharonian et al.l 120071 ). They are rather h igher than the values for the X-ray 
variability (e.g., 0.10-0.43 listed by lZhang et al.l (Il999l )). ( although much high er, and energy 
dependent, values of F^ar have been reported by MAGIC ( lAlbert et al.ll2007bl ) for Mrk 501). 



In Figure[7l our TeV light curve is plotted together with the Swift x-ray light curve (iFoschini et al 



20071 ). Unfortunately, Swift started observations two hours after CANGAROO stopped ob- 
serving each day, so there are no simultaneous observations. There is no evident correlation 
between the X-ray a nd the gamma-ray fluxes, with a time lag of approximately 9 hours. 
Foschini et al.l (120071 ) reported the X-ray flux increased by a factor of 5 in the 0.3-10 keV 
energy band without a large spectral change, and the highest peak of the Swift ligh t curve 
corresponds to the second large TeV flare observed by H.E.S.S. (IRaue et al.l l2006l ) . The 
X-ray flux in this outburst is not very high and it is not possible to detect time variation 
of the flux in the RXTE-ASM data. For observations in 2002 and 200 3, no correlation 



between the X-ray and the gamma-ray fluxes was reported by H.E.S.S. (lAharonian et al. 



2005bl ). however a strong corr e. 



ation during 2004 observations was presented in a prelim- 



inary analysis by iPunch et al.l (120071 ). On the other hand, Mrk 421 has displayed a good 



correlation between the X-ray and the gamma-ray fluxes, although recently no correlation 
results, including "Orphan flares", were reported. The simple one-zone SSC model is unable 
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to explain such complex behavior (e.g., iBlazejowski et al.l 120051). Some models ex plaining 



"Orphan flares" ha ve been proposed , for example, by iKusunose fc Takaharal (120061 ) for the 
leptonic model and iBottcher I (120051 ) for the hadronic model. 



During the outburst periods, the photon i ndex may change along with the flux variations 



as is the case for M rk 421 and Mrk 501 (e.g., lAharonian et al.l 120021 : iKrennrich et al.l 12002 



Albert et al.ll2007al jbll However, the average photon index we obtained here does not differ 



from the value that H.E.S.S. reported previously (lAharonian et al. 



2005a) within errors, 



and is co nsistent with those subs equently reported by H.E.S.S. (lAharonian et al.l 120071 ) and 
MAGIC ( iMazin &: Lindford 120071 ). The reason for this difference in the variation of photon 
index between PKS 2155—304 and Mrk 421 is still an open question. Although the error 
on the spectrum is enormous, a spectral index after correction for the absorption by the 
extra-galactic background light is calculated to be F = 1.6 ± 0.5 and F = 1.3 ±0.5 using the 



baseline model and the fast evolution model by lStecker &: ScuUyl (120061 ) . respectively. 



6. Conclusion 

We observed the southern HBL PKS 2155-304 from 2006 July 28 (MJD 53944) to 
August 2 with the CANGAROO-HI imaging Cherenkov telescopes. During the VHE high 
state we detected a signal at 4.8cr significance above 660 GeV in the total effective exposure 
time of 25.1 hours. The time averaged integral flux above 660 GeV is (1.6±0.3sf„j ±0.5syst) x 
10~^^ cm~^ sec~^ which corresponds to ~ 45% of the flux observed from the Crab nebula. The 
intra-night time variations of the flux were seen in our observations. Follow-up observations 
two weeks later indicated the source activity had decreased to lower than ~ 20% Crab flux.. 
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Table 1. Summary of observations for PKS 2155-304 from 2006 July 28 (MJD 53944) to 
August 2 (MJD 53944) and the follow-up observations from August 17 (MJD 53964) to 25 

(MJD 53972) ^ 







Begin 


End 


^obs 


z 


rtr 


Date 


Ntel 


[MJD] 


[MJD] 


[hrs] 


n 


[Hz] 


July 28 


2 


53944.514 


53944.663 


3.6 


29.8 


18.7 




3 


53944.668 


53944.830 


3.9 


20.4 


12.1 


July 29 


2 


53945.522 


53945.660 


3.1 


27.6 


14.9 




3 


53945.679 


53945.762 


2.0 


12.1 


6.2 


July 30 


2 


53946.535 


53946.656 


2.9 


26.1 


20.1 




3 


53946.668 


53946.833 


4.0 


22.2 


12.5 


July 31 


2 


53947.576 


53947.653 


1.8 


19.5 


19.8 




3 


53947.666 


53947.827 


3.9 


21.7 


11.6 


Aug. 2 


3 


53949.658 


53949.822 


3.9 


21.5 


11.9 


July 28-Aug.2 


2 






11.4 


26.6 


18.2 




3 






17.6 


20.4 


11.4 


Aug. 17-25 


3 






19.1 


20.9 


10.9 



^Observation date, the number of used telescopes, Nt^i, begin and 
end time of each observation in MJD, are summarized from the first 
column to the fourth column. In the following three columns, observa- 
tion time, tobs-^ average zenith angle, z, and average trigger rate, rtr are 
shown. Combined data between July 28 and August 2 are also shown. 
For the follow-up observations, only combined data are summarized. 
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Table 2. Summary of results for PKS 2155-304 from 2006 July 28 (MJD 53944) to 
August 2 (MJD 53944) and the follow-up observations from August 17 (MJD 53964) to 25 

(MJD 53972) ^ 



Date 






//,,[lirs] 






A' 




_F(>6G0G(A0 ^ 


July 28 


2 


12.4 


3.4 


675 


545 


130 ± 29 


4.6 


3.3 ±0.7 




3 


8.0 


3.5 


157 


132 


25± 14 


1.8 


1.0 ±0.6 


July 29^ 


2 


10.2 


2.2 


338 


340 


-2.4 ±20 


-0.12 


-0.10 ±0.8 




3 


4.1 


0.9 


39 


32 


7.1 ±6.1 


1.2 


1.3 ± 1.1 


July 30 


2 


13.2 


2.8 


434 


423 


11 ± 23 


0.48 


0.31 ± 1.0 




3 


8.1 


3.5 


235 


135 


100 ± 16 


6.1 


3.9 ±0.6 


July 31 


2 


13.0 


1.7 


229 


224 


4.8 ± 17 


0.29 


2.9 ± 1.0 




3 


7.6 


3.6 


156 


135 


21 ± 15 


1.4 


0.81 ± 0.57 


Aug. 2 


3 


7.8 


3.5 


142 


133 


9.4 ± 14 


0.66 


0.38 ± 0.57 


July 28-Aug.2 


2 


12.1 


10.1 


1676 


1532 


144 ± 56 


2.5 


1.3 ±0.5 




3 


7.4 


15.0 


729 


551 


178 ± 36 


5.0 


1.7 ±0.3 


Aug. 17-25 


3 


7.4 


17.1 


625 


637 


-12 ±36 


-0.34 


-0.12 ±0.34^1 



'^Following the observation date and the number of telescopes, average shower rate, Vgh, and live 
time, tiiy, are shown in the third and fourth columns, respectively. Nqn and Nqff in the fifth and 
sixth columns are the number of on-source and off-source events, respectively, where the off-source 
data from July 28 through August 2 were summed up for five nights and normalized to the live time. 
Seventh cohimn indicates the number of excess events derived by subtraction of off-source events from 
on-source events, and their significance are shown in the eighth column. Integral fiux above 660 GeV 
are calculated and are shown in the last column. Combined results from July 27 to August 2 and 
from August 17 and 25 are also presented, respectively. 

'^in unit of xl0~^^ cm^^ sec^^ 

'^affected by clouds 

^Threshold energy is 680 GeV. 
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Fig. 1. — Distributions of shower rate after image cleaning for 2fold{Mt panel) and 
3/oZd(riglit panel). Solid and hatched histograms show ones for zenith angles less than 
30° and larger than 30°, respectively. Prom our analysis we exclude the data whose shower 
rate are lower than 10 Hz, 8 Hz, and 5 Hz for 2fold-SZ, 2fold-LZ, and Sfold, respectively. 
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Fig. 2. — Distributions of squared angular distances 6*^ of 2fold and 3fold data for 
PKS 2155— 304 in the outburst period, obtained after the FD cut. (a):On-source data (error 
bars) and off-source data (hatched histogram) for 2fold are overlaid, where the latter is 
normalized to the live time (b): Closed circles show the excess events of the on-source above 
the off-source level for 2fold. The hatched histogram represents a PSF normahzed to the 
number of excess events in 9^ < 0.25. (c):Same as (a) for 3 fold. (d):Same as (b) for 3 fold. 
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Fig. 3. — Fisher Discriminant (FD) distribution for all events after 9'^ cut in the outburst 
period. Closed circles with error bars are on-source events {9^ < 0.06 for "if old data and 9"^ < 
0.14 for 2fold), and a solid histogram indicates off-source events. The closed squares show 
the distribution of on— off events which are gamma-ray candidates, and a dashed histogram 
is an expected distribution from Monte Carlo gamma-ray events normalized to the number 
of excess events. 
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Fig. 4. — Time-averaged differential energy spectrum of PKS 2155—304 between July 28 and 
August 2. Data plotted by closed circles are obtained from all the data at smaller zenith 
angle z < 30° {3fold-SZ and 2fold-SZ). The solid line represents the best fit to these data 
assuming a single power law spectrum. For comparison, the time av erage differential energy 
spectra in early July 28 in 200 7 reported by lAharonian et al.l (120071 ) during a giant flare and 
in 2002 and 2003 reported by lAharonian et al.l (l2005al ) at a low state are drawn by dashed 
and dotted lines, respectively. 
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Fig. 5. — Daily average integral flux of PKS 2155— 304 using all the data above 660 GeV 
between 2006 July 28 and August 2. Dashed line indicates an averaged integral flux during 
this observation period (~ 45 % Crab flux level). 
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Fig. 6. — Light curve of PKS 2155— 304 using all the data at zenith angle less than 30° 
between 2006 July 28 and August 2, expressed by the integral flux above 660 GeV. Closed 
triangles and closed circles indicate the results from 2fold-SZ and Sfold-SZ datasets, 
respectively. Dashed line indicates an average integral flux during this observation period. 
The bin width is 40 minutes. The shaded areas indicate the H.E.S.S. observation periods. 
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Fig. 7. — Comparison of the light curves between VHE gamma-ray and X-ray bands, (a) A 
nightly average integral flux of PKS 2155-304 above 660 GeV obtained by the CANGAROO- 
III given in unit of cm~^s~^, which is the same as in Fig. [51 An upper limit of the integral 
flux between August 17 and 25 is also plotted, which is an extrapolate value from the actual 
upper limi t above 680 GeV. (b) X-ray [0.3-10 keV] counts in unit of count s~^ from the XRT 
on Swift (IFoschini et al.l 120071 ). 



